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Abstract 

We have reconsidered the atmospheric neutrino anomaly in light of the 1289-day data from 
Super-Kamiokande contained events and from Super-Kamiokande and MACRO up-going muons. 
We have reanalysed the proposed solution to the atmospheric neutrino anomaly in terms of non- 
standard neutrino-matter interactions (NSI) as well as the standard — > Vr oscillations (OSC). 
Our statistical analysis shows that a pure NSI mechanism is now ruled out at 99%, while the stan- 
dard Vr OSC mechanism provides a quite remarkably good description of the anomaly. We 
therefore study an extended mechanism of neutrino propagation which combines both oscillation 
and non-standard neutrino-matter interactions, in order to derive limits on flavour-changing (FC) 
and non-universal (NU) neutrino interactions. We obtain that the off-diagonal flavour-changing 
neutrino parameter e and the diagonal non-universality neutrino parameter e' are confined to 
—0.05 < e < 0.04 and \e'\ < 0.17 at 99% CL. These limits are model independent and they are 
obtained from pure neutrino-physics processes. The stability of the neutrino oscillation solution 
to the atmospheric neutrino anomaly against the presence of non-standard neutrino interactions 
establishes the robustness of the near-maximal atmospheric mixing and massive-neutrino hypoth- 
esis. The best agreement with the data is obtained for /S.m? = 2.4 x 10~^ eV^, siv?{26) = 0.99, 
e = —9.1 X 10~^ and e' = —1.9 x 10~^, although the function is quite flat in the e and e' directions 
for e, e' — > 0. A revised analysis which takes into account the new 1489-day Super-Kamiokande 
and final MACRO data is presented in the appendix; the determination of Am^ and 6 is essentially 
unaffected by the inclusion of the new data, while the bounds on e and e' are strongly improved 
to -0.03 <e< 0.02 and le'l < 0.05 at 99.73% CL. 
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I. INTRODUCTION 

The experimental data on atmospheric neutrinos show, in the muon-type 

events, a clear deficit which cannot be accounted for without invoking non-standard neutrino 
physic. This result, together with the solar neutrino anomaly g is very important since 
it constitutes a clear evidence for physics beyond the Standard Model. Altogether, the 
simplest joint explanation for both solar and atmospheric anomalies is the hypothesis of 
three-neutrino oscillations 0. 

There are however many attempts to account for neutrino anomalies without oscilla- 
tions [2i]- Indeed, in addition to the simplest oscillation interpretation 01, the solar 
neutrino problem admits very good alternative explanations, for example based on tran- 
sition magnetic moments or non-standard neutrino interactions (NSI) Likewise, 
several such alternative mechanisms have been postulated to account for the atmospheric 
neutrino data such as the NSI [12] or the neutrino decay hypotheses 



In contrast to the solar case, the atmosp^ 



leric neutrino anomaly is so well reproduced by 



the ^ osculation hypothecs (OSC) Q that one can use the robus;ness of this 

interpretation so as to place stringent limits on a number of alternative mechanisms. 

Among the various proposed alternative interpretations, one possibility is that the neutri- 
nos posses non-standard interactions with matter, which were shown to provide a good de- 
scription of the contained event data sample [l^. Such non-standard interactions [3,E3,Ell 
can be either flavour-changing (FQ or non-universal ( NU) , and arise naturally in theoreti- 
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cal models for massive neutrinos 22 , |23|, |2J, |25|, |26|, |27|, |28[ . This mechanism does not even 
require a mass for neutrinos |2 
in most models 



although neutrino masses are expected to be present 
251 . |26L l27l . I28l lio^ . It is therefore interesting to check whether the 
atmospheric neutrino anomaly could be ascribed, completely or partially, to non-standard 
neutrino-matter interactions. In Refs. 0, 1^ the atmospheric neutrino data have been 
analysed in terms of a pure z/^ conversion in matter due to NSI. The disappearance 

of V/j, from the atmospheric neutrino flux is due to interactions with matter which change 
the flavour of neutrinos. A complete analysis of the 52 kton-yr Super-Kamiokande data was 
given in Ref. [s^l- It included both the low-energy contained events as well as the higher 
energy stopping and through-going muon events, and showed that the NSI solution was 
acceptable, although the statistical relevance was low. Compatibility between the data and 
the NSI hypothesis was found to be 9.5% for relatively large values of flavour-changing and 
non-universality parameters^ . 

In the present paper we will use the latest higher statistics data from Super-Kamiokande 
(79 kton-yr) [2] and MACRO data in order to briefly re-analyse the atmospheric data 



^ For more exotic attempts to explain the neutrino anomalies see 0, 0, Q| . 

^ For another analysis showing low confidence for a dominant NSI in atmospheric neutrinos, see [32| . 
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within the oscillation hypothesis. We show that the oscillation description has a high signif- 
icance, at the level of 99% for the Super-Kamiokande data, and of 95% when the MACRO 
through-going muons data are also added to the analysis. We then show that the new data 
rule out the NSI mechanism as the dominant conversion mechanism. The goodness of the 
fit (GOF) is now lowered to 1%. This clearly indicates that a pure NSI mechanism can not 
account for the atmospheric neutrino anomaly. 

However, the possibility that neutrinos both posses a mass and non-standard interactions 
is an intriguing possibility. For example in models where neutrinos acquire a mass in see-saw 
type schemes the neutrino masses naturally come together with some non-diagonality of the 
neutrino states ji^. Alternatively, in supersymmetric models with breaking of R parity [3| 
neutrino masses and flavour-changing interactions co-exist^. This in turn can induce some 
amount of flavour-changing interactions. The combined mechanism of oscillations (OSC) 
together with NSI may be active in depleting the atmospheric z/^ flux, and therefore it can 
provide an alternative explanation of the deficit. Since the atmospheric neutrino anomaly 
is explained remarkably well by z/^ — > Ur oscillations, while pure NSI cannot account for the 
anomaly, this already indicates that NSI can be present only as a sub-dominant channel. The 
atmospheric neutrino data can therefore be used as a tool to set limits to the amount of NSI 
for neutrinos. These limits are obtained from pure neutrino-physics processes and are model 
independent, since they do not rely on any specific assumption on neutrino interactions. 
In particular they do not rely on any SU{2)l assumption relating the flavour-changing 
neutrino scattering off quarks (or electrons) to interactions which might induce anomalous 
tau decays [ssl or suffer from QCD uncertainties. In the following we will show that, from 
the analysis of the full set of the latest 79 kton-yr Super-Kamiokande and the MACRO 
data on up-going muons atmospheric neutrino data, FC and non-universal neutrino 
interactions are constrained to be smaller than 5% and 17% of the standard weak neutrino 
interaction, respectively, without any extra assumption. 

The plan of the paper is the following. In Sec. |n] we briefly describe the theoretical 
origin of neutrino NSI in Earth matter. In Sec. II I II we briefly summarize our analysis of the 
atmospheric neutrino data in terms of z/^ vacuum oscillations. In Sec. IIVI we update 

our analysis for the pure NSI mechanism, and we show that the latest data are able to rule it 
out as the dominant — > z/,- conversion mechanism for atmospheric neutrinos. In Sec. IVlwe 
therefore investigate the combined situation, where massive neutrinos not only oscillate but 
may also experience NSI with matter. In this section we derive limits to the NSI parameters 
from the atmospheric neutrino data. In Sec. ED we present our conclusions. 



The NSI may, however, be rather small 



H- 
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II. THEORY 



Generically models of neutrino mass may lead to both oscillations and neutrino NSI in 
matter. Here we sketch two simple possibilities. 



A. NSI from neutrino-mixing 

The most straightforward case is when neutrino masses follow from the admixture of 
i=o=inglet neutral heavy lepton= as, tor example, in seesaw schemes 0. These contain 

SU{2) ® U{1) singlets with a gauge invariant Majorana mass term of the type MRijV^v^ 
which breaks total lepton number symmetry. The masses of the light neutrinos are obtained 
by diagonalizing the mass matrix 

in the basis z/, z/*^, where D is the standard SU{2) ® U{1) breaking Dirac mass term, and 



(1) 



221. 



Mji = is the large isosinglet Majorana mass and the Mlvv term is an iso-triplet 
In S'O(IO) models the first may arise from a 126 vacuum expectation value, while the latter 
is generally suppressed by the left-right breaking scale. Ml oc 1/Mr. 

In such models the structure of the associated weak currents is rather complex [s^l- The 
first point to notice is that the isosinglets, presumably heavy, will mix with the ordinary 
isodoublet neutrinos in the charged current weak interaction. As a result, the mixing matrix 
describing the charged leptonic weak interaction is a rectangular matrix K which may 
be decomposed as 

K={Kl.Kh) (2) 

where and Kh are 3x3 matrices. The corresponding neutral weak interactions are 
described by a non-trivial matrix 



22 



P = K^K. (3) 

In such models non-standard interactions of neutrinos with matter are of gauge origin, 
induced by the non-trivial structures of the weak currents. Note, however, that since the 
smallness of neutrino mass is due to the seesaw mechanism M^eff = Mi — DM^^D^ the 
condition 

Ml < Mr (4) 

the magnitude of neutrinos NSIs is expected to be neghgible. 

However the number m of SU{2) ® f/(l) singlets is completely arbitrary, so that one may 
consider the phenomenological consequences of models with Majorana neutrinos based on 
any value of m. In this case one has 3(1 -|- m) mixing angles 6ij and the same number of 
CP violating phases (pij characterizing the neutrino mixing matrix K 22, "s^. This number 
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far exceeds the corresponding number of parameters describing the charged current weak 
interaction of quarks. The reasons are that (i) neutrinos are Majorana particles so that their 
mass terms are not invariant under rephasings, and (ii) the isodoublet neutrinos mix with 
the isosinglets. For m < 3, 3 — m neutrinos will remain massless, while 2m neutrinos will 
acquire Majorana masses but may have non-zero NSI. For example, in a model with m = 1 
one has one light neutrino and one heavy Majorana neutrino in addition to two massless 
neutrinos whose degeneracy is lifted by radiative corrections. 

In contrast, the case m > 3 may also be interesting because it allows for an elegant 
way to generate neutrino masses without a superheavy scale, such as in the seesaw case. 
This allows one to enhance the allowed magnitude of neutrino NSI strengths by avoiding 
constraints related to neutrino masses. As an example consider the following extension of 
the lepton sector of the SU(2) ® U (1) theory: let us add a set of two 2-component isosinglet 
neutral fermions, denoted z/^j and Si, in each generation. In this case one can consider the 
9x9 mass matrix [2^ 

"0 DO 

M 

(in the basis z/, z/'^, S). The Majorana masses for the neutrinos are determined from 



(5) 



(6) 



In the limit /i — > the exact lepton number symmetry is recovered and will keep neutrinos 
strictly massless to all orders in perturbation theory, as in the Standard Model [2^. The 
propagation of the light (massless when /i — > 0) neutrinos is effectively described by an 
effective truncated mixing matrix Kr which is not unitary. This may lead to oscillation 
effects in supernovae matter, even if neutrinos were massless The strength 

of NSI is therefore unrestricted by the magnitude of neutrino masses, only by universality 
limits, and may be large, at the few per cent level. The phenomenological implications of 
these models have been widely investigated 



B. NSI from new scalar interactions 

An alternative and elegant way to induce neutrino NSI is in the context of unified super- 
symmetric models as a result of sup ersy mmetric scalar lepton non-diagonal vertices induced 
by renormalization group evolution 00. In the case of SU (5) the NSI may exist without 
neutrino mass. In 50(10) neutrino masses co-exist with neutrino NSI. 

An alternative way to induce neutrino NSI without invoking physics at very large mass 
scales is in the context of some radiative models of neutrino masses In such models 

NSI may arise from scalar interactions. 
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FIG. 1: Diagram generating neutrino mass in supersymmetry with explicitly broken R-parity. It 
illustrates the co-existence of OSC and NSI mechanisms used in Eq. (^Hl 

Here we focus on a more straightforward way to induce NSI based on the most general 
form of low-energy supersymmetry. In such models no fundamental principle precludes 
the possibility to violate R parity conservation j^l explicitly by renormalizable (and hence 
a priori unsuppressed) operators such as the following extra L violating couplings in the 
superpotential 

\jkLiLjE1 (7) 
K,kUQ,Dl (8) 

where L, Q, and D*^ are (chiral) superfields which contain the usual lepton and quark 
SU{2) doublets and singlets, respectively, and i^j^k are generation indices. The couplings 
in Eq. (|7j) give rise at low energy to the following four-fermion effective Lagrangian for 
neutrinos interactions with (i-quark including 

Leff = -2\/2GF^iapVLal^VLpdRl^dR a,p = e,n,T, (9) 

a,f3 

where the parameters ^ai3 represent the strength of the effective interactions normalized to 
the Fermi constant Gp- One can identify explicitly, for example, the following non-standard 
flavour-conserving NSI couplings 

£ =V (11) 



and the FC coupling 



f - V ^'^^^^'^^^ (12) 

J ^ ^ IjL 



where rriq.^ are the masses of the exchanged squarks and j = 1,2,3 denotes dL,SL,bL, 
respectively. Likewise, one can identify the corresponding flavour-changing NSI. The ex- 
istence of effective neutral current interactions contributing to the neutrino scattering off 



-2 




0.6 0.7 0.8 0.9 1 

sin^29 

FIG. 2: Allowed regions in the IS.m?-s\v?{29) parameter space for the pure Vr oscillation 

mechanism. The shaded areas refer to the 90%, 95%, 99% and 99.73% CL with 2 parameters. 
The best fit point is indicated by a star. Both Super-Kamiokande and MACRO data have been 
included. 

d quarks in matter, provides new flavour-conserving as well as flavour-changing terms for 
the matter potentials of neutrinos. Such NSI are directly relevant for atmospheric neutrino 
propagation. As a final remark we note that such neutrino NSI are accompanied by non-zero 
neutrino masses, for example, induced by loops such as that in Fig. [H The latter lead to vac- 
uum oscillation (OSC) of atmospheric neutrinos. The relative importance of NSI and OSC is 
model-dependent. In what follows we will investigate the relative importance of NSI-induced 
and neutrino mass oscillation induced (OSC-induced) conversion of atmospheric neutrinos 
allowed by the present high statistics data. 



III. VACUUM OSCILLATION HYPOTHESIS 

We first briefly report our updated results for the usual —>■ v-r vacuum oscillation 
channel. For definiteness we confine to the simplest case of two neutrinos, in which case CP 
is conserved in standard oscillations'^. The evolution of neutrinos from the production point 
in the atmosphere up to the detector is described by the evolution equation: 

(13) 



:: 













In L-violating oscillations there is in principle CP violation due to Majorana phases. 
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FIG. 3: Zenith-angle distributions for the Super-Kamiokande and MACRO data sets, together 
with our predictions in the absence of oscillation (thick solid line) and the predictions for the best 
fit points for each data set in the different Vj- transition channels: pure oscillation (thin solid 

line), pure NSI (dashed line) and the hybrid oscillation + NSI mechanism (dot-dashed line). The 
errors displayed in the experimental points are statistical only. 



where the Hamiltonian which governs the neutrino propagation can be written as: 



H 



Hut Ht- 




(14) 



In Eq. fll4|) Am^ is the squared-mass difference between the two neutrino mass eigenstates 
and the rotation matrix Re is simply given in terms of the mixing angle 9 by 



cos 9 sin 6' 
- sin 9 cos 9 



(15) 



The oscillation probability for a neutrino which travels a path of length L is therefore: 

Am^L 



sin^(2^) sin^ 1.27 



(16) 



where Am^, L and Ey are measured in eV^, Km and GeV, respectively. 

The calculation of the event rates and the statistical analysis is performed according to 
Ref. 01 . In the present analysis we include the full set of 79 kton-yr Super-Kamiokande 



8 







Vr oscillations 


NSI hypothesis 


Data Set 


d.o.f. 


Am^ [eV2] sin2(2e) Xosc GOF 


e e' xlcGOF 


SK Sub-GeV 
SK Multi-GeV 
SK Stop-^ 
SK Thru-/x 
MACRO 


10-2 
10-2 
5-2 
10-2 
10-2 


2.2 X 10-3 1.00 4.1 84% 
2.1 X 10-3 0.94 4.2 84% 
3.0 X 10-3 0.99 0.7 88% 

6.3 X 10-3 0.78 5.3 73% 
1.3 X 10-3 1.00 11.6 17% 


0.196 0.010 5.1 75% 
0.667 0.431 4.2 84% 
0.697 0.317 2.5 48% 
0.041 0.138 5.7 68% 
0.020 0.046 6.6 58% 


SK Contained 
SK Upgoing 
SK Cont+Stop 
Thru-;U 


20-2 
15-2 
25-2 
20-2 


2.1 X 10-3 1.00 8.8 96% 

3.2 X 10-3 0.94 6.5 92% 
2.5 X 10-3 0.99 10.0 99% 
3.0 X 10-3 0.95 18.1 45% 


0.667 0.138 10.9 90% 
0.041 0.144 16.5 22% 
0.697 0.331 15.3 88% 
0.018 0.058 21.1 28% 


SK 

SK+MACRO 


35-2 
45-2 


2.7 X 10-3 0.97 16.2 99% 
2.5 X 10-3 0.96 28.7 95% 


0.536 0.611 53.1 1% 
0.513 0.667 67.6 1% 



TABLE I: Minimum values and best-fit points for the various atmospheric neutrino data sets 
considered in the analysis and for two different neutrino conversion mechanisms: pure — > Vt 
vacuum oscillation (OSC) and pure non-standard neutrino-matter interactions (NSI). 



data and the latest MACRO data on upgoing muons [2^. The results of the fits are 
shown in TableHl the best fit point is Am^ = 2.7 x 10-3 eV^ and sin^ 29 = 0.97 with a GOF 
of 99% when only Super-Kamiokande data are considered. The inclusion of MACRO lowers 
slightly the GOF to 95% but practically does not move the best fit point, which in this case 
is Am2 = 2.5 x 10-3 eV^ and sin^ 26 = 0.96. 

Fig.|21shows the allowed region in the plane (sin^ 29, Am^), and Fig.Elreports the angular 
distributions of the Super-Kamiokande data sets and the same distributions calculated for 
the best fit point. The agreement between the data and the calculated rates in presence of 
oscillation is remarkable, for each data sample. The same occurs also for the MACRO data 
set. 

From this analysis we can conclude that the z/^ —>■ v-r oscillation hypothesis represents a 
remarkably good explanation of the atmospheric neutrino anomaly (see also Refs. 17. isj]). 



IV. NON STANDARD NEUTRINO INTERACTIONS 



Let us re-analyze the interpretation of the atmospheric neutrino anomaly in terms of pure 



non-standard interactions of neutrinos with matter 
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30 



3l|. In this case, neutrinos are 



assumed to be massless and the v., 



Vt conversion is due to some NSI with the matter 



which composes the mantle and the core of the Earth. The evolution Hamiltonian can be 
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FIG. 4: Function a of Eq. (|23|) and the relevant product (aL) which enters in the pure NSI 
transition probabihty of Eq. (|2'2j) . plotted as a function of the cosine of the Earth's zenith angle rj. 



written as 
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ii = ±y/2GFNf{r 



e,. 



(17) 



where the sign + (— ) holds for neutrinos (antineutrinos) and and e'^^ parametrize the 
deviation from standard neutrino interactions: \/2GpNf{r)ei, is the forward scattering am- 
plitude of the FC process + f — > i^t + f and \/2 G pN f{r)e'jj represents the difference 
between the Vr + f and the + f elastic forward scattering amplitudes. The quantity 
Nf{r) is the number density of the fermion / along the path r of the neutrinos propagating 
in the Earth. To conform to the analyses of Ref. 0], we set our normalization on these 
parameters by considering that the relevant neutrino interaction in the Earth occurs only 
with down-type quarks. 

In general, an equation analogous to Eq. (fTTjl holds for anti-neutrinos, with parameters 
Ep and e'p. For the sake of simplicity, we will assume here and in the following = Ep = e 
and el = e'p = e' . It is therefore useful to introduce the following variables (F, y?) instead of 

e = Fsm{2(p), 



-=Fcos(2<^), 



(18) 



or, equivalently: 



£'74 + £2 



■ arctan 



(19) 



2 V^V2. 

With the use of the variables F and 6, the evolution Hamiltonian Eq. ()17|1 can be cast in a 
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FIG. 5: Allowed regions in the e-e' parameter space for the pure v-r NSI mechanism and for 

different sets of experimental data. The shaded areas refer to the 90%, 95%, 99% and 99.73% CL 
with 2 parameters. For each panel, the best fit point is indicated by a star. 



form which is analogous to the standard oscillation one: 

n = ±V2GFNf{r)Fn^(^^^ ^^Rt,, (20) 
where R<^ assumes the structure of a usual rotation matrix with angle y?: 

R,= f ""^V (21) 

\ — sin (y9 COS ip I 

The transition probabilities of u^^ Vr (i^fj. —>■ ^t) are obtained by integrating Eq. ()2()|1 
along the neutrino trajectory inside the Earth. For the Earth's density profile we employ 
the distribution given in and a realistic chemical composition with proton/nucleon ratio 
0.497 in the mantle and 0.468 in the core Although the integration is performed 

numerically, the transition probability can be written exactly in a simple analytical form as 

Pu.^u. = Pu.^p. = sin^ {2ip) sin^ {aF L), (22) 
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where 

a = V2Gf {Nf) (23) 

and (Nf) is the mean value of Nf{r) along the neutrino path. Note that the analytical form 
in Eq. holds exactly despite the fact that the number density Nf{r) varies along the 
path. The quantity a and the relevant product aL which enters the transition probability 
in Eq. fl22j) are plotted in Fig. |3] as a function of the zenith angle t] and calculated for the 
Earth's profile quoted above. From Fig.|3]it is clear the sharp change from the mantle to the 
core densities which occurs for cost] ~ 0.84. Notice that the transition probability P^^^i,^ 
(Pp^^p^) is formally the same as the expression for vacuum oscillation Eq. ()16|) with the angle 
ip playing a role of mixing angle analogous to the angle 9 for vacuum oscillations. In the 
other hand, in the factor which depends on the neutrino path L, the parameter F formally 
replaces Am^. However, in contrast to the oscillation case, there is no energy dependence 
in the case of NSI IH H E] . 

The result of the fits to the Super-Kamiokande and MACRO data are reported in Fig. 
and again in Table HI As already discussed in Ref. the NSI mechanism properly accounts 
for each Super-Kamiokande data set separately, as well as the MACRO upgoing muons 
data. Moreover it succeeds in reconciling together the sub-GeV, multi-GeV and stopping 
muons data sets. However, the NSI cannot account at the same time also for the through- 
going muons events, mainly because the NSI mechanism provides an energy independent 
conversion probability, while the upgoing muon events, which are originated by higher energy 
neutrinos, require a suppression which is smaller than the one required by the other data 



sets 



12, 



30, 



3l|. This effect is clearly visible in two ways. First, from Fig. where we can 
see that the allowed regions for SK contained + stopping-/^ events (upper-right panel) and 
for SK + MACRO through-going /z events (lower-left panel) are completely disjoint even 
at the 99.73% CL. In addition, from the angular distribution of the rates shown in Fig. 01 
where the angular distribution for upgoing muons calculated for the best fit point of the 
pure NSI mechanism clearly shows too a strong suppression, especially for horizontal events. 
The global analysis of Super-Kamiokande and MACRO data has a very low GOF, only 1%: 
this now allows us to rule out at 99% the pure NSI mechanism as a possible explanation of 
the atmospheric neutrino anomaly. 



V. COMBINING THE OSC AND NSI MECHANISMS 

Let us now consider the possibility that neutrinos are massive and moreover posses non- 
standard interactions with matter. As mentioned in Sec. m this may be regarded as generic 
in a large class of theoretical models. In this case, their propagation inside the Earth is 
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sin^20 



FIG. 6: Allowed regions in the Am^-siia^ {29) parameter space for the hybrid OSC + NSI mecha- 
nism. In each panel, the value of the NSI parameter F is fixed, while the other NSI parameter ip 
is integrated out. The last panel shows the allowed region when both F and (p are integrated out. 
The shaded areas refer to the 90%, 95%, 99% and 99.73% CL with 3 parameters, and the best fit 
point is indicated by a star. Both Super-Kamiokande and MACRO data have been included. 



governed by the following Hamiltonian 



H 



-1 
1 



Rl±V2GFNf(r)FR, 



-1 
1 



(24) 



where and R;^ are the mixing matrices defined in Eqs. (|TH|) and respectively. The 
NSI term in the Hamiltonian has an effect which is analogous to the presence of the effective 
potentials for the propagation in matter of massive neutrinos, a situation which leads to 
the MSW oscillation mechanism 47[. Also in the case of Eq. ()24|) neutrinos can experience 



matter-induced oscillations, due to the fact that z/„'s and i/^'s can have both flavour-changing 
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FIG. 7: Allowed regions in the F-ip parameter space for the hybrid OSC + NSI mechanism. The 
parameters Arn^ and sin^(20) are integrated out. The shaded areas refer to the 90%, 95%, 99% 
and 99.73% CL with 2 parameters, and the best fit point is indicated by a star. The top panel 
shows the behaviour of the & function of the NSI parameter F when cp is also integrated out. 

Both Super-Kamiokande and MACRO data have been included. 



and non-universal interaction with the Earth matter. 

Since the Earth's matter profile function Nf{r) is not constant along the neutrino propa- 
gation trajectories, the Hamiltonian matrices calculated at different points inside the Earth 
do not commute. This leads to a non trivial evolution for the neutrinos in the Earth and a 
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8 

FIG. 8: Allowed regions in the e-e' parameter space for the hybrid OSC + NSI mechanism. The 
parameters Am^ and sin^(20) are integrated out. Notice that both positive and negative values 
of e and e' are shown. The shaded areas refer to the 90%, 95%, 99% and 99.73% CL with 2 
parameters, and the best fit point is indicated by a star. Both Super-Kamiokande and MACRO 
data have been included. 

numerical integration of the Eq. (fT!^ with the Hamiltonian of Eq. is needed in order to 
calculate the neutrino and anti- neutrino transition probabilities P^^^ur and Pp^^p^. 

The transition mechanism depends on four independent parameters: the neutrino 
squared-mass difference Am^, the neutrino mixing angle 6, the PC parameter e and the 
NU parameter e' (or, alternatively, the F and ip parameters for the NSI sector). In our 
analysis we will use the F and if parameters, which prove to be more useful, and then ex- 
press the results, which we will obtain for these two parameters, in terms of the e and e' 
parameters, which have a more physical meaning. 

As a first step, we can use the symmetries of the Hamiltonian in order to properly define 
the intervals of variation of the parameters. Since H in Eq. ()24|1 is real and symmetric, the 
transition probabilities are invariant under the following transformations: 
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FIG. 9: Behaviour of the as a function of the flavour-changing parameter e (top panel) and of 
the non-universal neutrino-interactions parameter e' (bottom panel), for the hybrid OSC + NSI 
mechanism. In each panel both the oscillation parameters (sin^(20) and Am^) and the undisplayed 
NSI parameter (e' in the top panel and e in the bottom one) are integrated out. Notice that both 
positive and negative values of e and e' are shown. 



• if ^ if + 11, 

• Am"^ -Am'^ and 6 6 + 7^/2, 

• F —>■ —F and ip^ip + TT/2. 

Under any of the above transformations the Hamiltonian remains invariant. Moreover even 
if the overall sign of the Hamiltonian changes this will have no effect in the calculation of 

Py^^^ur and Pp^^p^: 

• 9^9 + n/2 and Lp^Lp + 7r/2 (or: e ^ —e and e' — > — e'). 

Finally, if the sign of the non-diagonal entries in the Hamiltonian changes, again there is no 
effect in the neutrino/anti-neutrino conversion probabilities: 

• 9 ^ —9 and ip —ip (or: e — £:). 
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The above set of invariance transformations allows us to define the ranges of variation of 
the four parameters as follows: 



(a) O<0<7r/4, 

(b) < < TT , 

^ ' - ^ - ' (25) 

(c) Am^ > , 

(d) F > . 

Notice that, in contrast with the MSW mechanism, it is possible here, without loss of 
generality, to constrain both the mixing angle 6 inside the [0, 7r/4] interval keeping Am^ 



positive. There is no "dark side" |4^ in the parameter space for this mechanism^. In our 
analysis we will adopt the set of conditions of Eq. , implying that the neutrino squared- 
mass difference and mixing angle are confined to the same intervals as in the standard 
i^ti —>■ oscillation case, while the NSI parameters e and e' can assume independently both 
positive and negative values. We will actually find that the best fit point occurs for negative 
e and e' . 

Let us turn now to the analysis of the data and the presentation of the results. Here we 
perform a global fit of the Super-Kamiokande data sets and of the MACRO upgoing muon 
flux data in terms of the four parameters of the present combined OSC + NSI mechanism. 
As we have already seen in the previous sections, pure oscillation provides a remarkably good 
fit to the data, while the pure NSI mechanism is not able to reconcile the anomaly observed 
in the upgoing muon sample with that seen in the contained event sample. This already 
indicates that, when combining the two mechanisms of — > v-r transition, the oscillation 
will play the role of leading mechanism, while the NSI could be present at a subdominant 
level. 

As a first result, we quote the best fit solution: Am? = 2.4 x 10"^ eV^, sin^(26') = 0.99, 
e = —9.1 X 10~^ and e' = —1.9 x 10"^. The goodness of the fit is 94% (45 — 4 degrees 
of freedom). For the Am^ and sin^(2^) parameters, the best fit is very close to the best 
fit solution for pure oscillation (see Table |ll). This is a first indication that the oscillation 
mechanism is stable under the perturbation introduced by the additional NSI mechanism. 
It is interesting to observe that a small amount of FC could be present, at the level of less 
than a percent, while z/^ and z/^ interactions are likely to be universal. Moreover, the 
function is quite fiat in the e and e' directions for e, e' —* 0. 

We also display the effect of the NSI mechanism on the determination of the oscillation 
parameters by showing the result of the analysis in the Am^ and sin^(20) plane, for fixed 
values of the NSI parameters. 



^ We also notice that one can replace conditions (a) and (6) in Eq. (|25|l by (a') < 9 < tt and {b') < (f < 
7r/4. This implies that both e and e' are positive in this case. 
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Fig. ini shows the dependence of the allowed region in the Am^ and sin^(2^) plane for 
fixed values of the NSI parameters, in particular for fixed values of F irrespective of the 
value of <y9, which is "integrated out". Note that for F < 0.02 the allowed region is almost 
unaffected by the presence of NSI. For larger values the quality of the fit gets rapidly worse, 
however the position of the best fit point in the plane (sin^(26'), Am^) remains extremely 
stable. For F > 0.1 the 99% CL allowed region finally disappears. The last panel of Fig. IHl 
shows the allowed region when both F and if are integrated out. The region obtained is in 
agreement with the one obtained for pure oscillation case. We can therefore conclude that 
the determination of the oscillation parameters Am^ and sin^(26') is very stable under the 
effect of non-standard neutrino-matter interactions. 

We can now look at the results from the point of view of the NSI parameters. This will 
allow us to set bounds on the maximum allowed level of neutrino NSI. Fig. [7| shows the 
behaviour of the as a function of the F parameter, and the allowed region in the F and 
(f parameter space with Am^ and sin^(2^) integrated out. From the lower panel we see that 
the F parameter is constrained by the data to values smaller than ~ 0.09 at 99% CL, while 
the quantity ip is not constrained to any specific interval. When ip is also integrated out 
(upper panel of Fig. Ej) the number of free parameters is reduced to 1, and the upper bound 
on F improves to ~ 0.05. 

Looking at Fig. [7| and taking into account the definition of F and (p in terms of e and e' 
given in Eq. (fT^ . we see that the data constrain the maximum amount of FC and NU which 
is allowed (from F), but they do not fix their relative amount (through ^p). This information 
can be conveniently translated in the e and e' plane, as we show in Fig. |H| at 99% CL, the 
flavour-changing parameter e is confined to —0.05 < e < 0.04, while the non-universality 
parameter is bounded to |£:'| < 0.17. These are the strongest bounds which can be imposed 
simultaneously on both FC and NU neutrino-matter interactions, but it is also interesting to 
look at the separate behaviour of the with respect to either FC or NU-type neutrino NSI 
when the other type of interaction is also integrated out. This is illustrated in Fig. IHl where 
we see that the bounds on e and e' - now calculated with only 1 degree of freedom - are 
improved to —0.03 < e < 0.02 and \e'\ < 0.07. We also notice that the function is more 
shallow for e' than for e, indicating that the bound on FC interactions is more stringent 
than the one on NU interactions. 

This is the main result of our analysis, since it provides limits to non-standard neutrino 
interactions which are truly model independent, since they are obtained from pure neutrino- 
physics processes. In particular they do not rely on any relation between neutrinos and 
charged lepton interactions. Therefore our bounds are totally complementary to what may 
be derived on the basis of conventional accelerator experiments jl^. Note that although the 
above bounds of neutrino-matter NSI were obtained simply on the basis of the quality of 
present atmospheric data, they are almost comparable in sensitivity to the capabilities of a 
future neutrino factory based on intense neutrino beams from a muon storage ring [Is^ . 
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VI. CONCLUSIONS 



In this paper we have analysed the most recent and large statistic data on atmospheric 
neutrinos (Super-Kamiokande and MACRO) in terms of three different mechanisms: (i) pure 
OSC Vn Vj. oscillation; (ii) pure NSI transition due to non-standard neutrino- 

matter interactions (flavour-changing and non-universal); (iii) hybrid OSC + NSI i/r 
transition induced by the presence of both oscillation and non-standard interactions. 

The pure oscillation CclSG, clS IS well known, provides a remarkably good fit to the ex- 
perimental data, and it can be considered the best and most natural explanation of the 
atmospheric neutrino anomaly. In this updated analysis, we obtain the best fit sohition for 
Am^ = 2.5 X 10-3 eV^ and sin^ 29 = 0.96, with a goodness-of-fit of 95% (Super-Kamiokande 
and MACRO combined). 

In contrast, the pure NSI mechanism, mainly due to its lack of energy dependence in the 
transition probability, is not able to reproduce the measured rates and angular distributions 
of the full data sample because it spans about three orders of magnitude in energy. The 
data clearly show the presence of an up-down asymmetry and some energy dependence. 
With the increased statistics of the data presently available it is now possible to rule out 
this mechanism at 99% as a possible explanation of the atmospheric neutrino data. 

We have therefore investigated a more general situation: the possibility that massive 
neutrinos also possess some amount of flavour-changing interactions with matter, as well 
as some difference in the interactions between i/^'s and i/^-'s. The global analysis of the 
Super-Kamiokande and MACRO data shows that the oscillation hypothesis is very stable 
against the possible additional presence of such non-standard neutrino interactions. The 
best fit point is obtained for Am^ = 2.4 x lO'^ eV^ sin^(2^) = 0.99, e = -9.1 ■ lO'^ 
and e' = — 1.9 x 10^'^ with a goodncss-of-fit of 94% (45 — 4 degrees of freedom). A small 
amount of FC could therefore be present, at the level of less than a percent, while z/^ and 
i/r interactions are likely to be universal. In addition the function is rather fiat in the e 
and e' directions for e,e' ^ and NSI can be tolerated as long as their effect in atmospheric 
neutrino propagation is subdominant. 

From the analysis we have therefore derived bounds on the amount of flavour-changing 
and non-universahty allowed in neutrino-matter interactions. At the 99% CL, the flavour- 
changing parameter e and the non-universality parameter e' are simultaneously confined to 
—0.05 < £ < 0.04 and |£:'| < 0.17. The bounds on fiavour-changing interactions is stronger 
than the one which applies on universality violating ones. These bounds on non-standard 
neutrino interactions do not rely on any assumption on the underlying particle physics 
model, as they are obtained from pure neutrino-physics processes. They could be somewhat 
improved at a future neutrino factory based on intense neutrino beams from a muon storage 
ring. 

Note in particular that the bounds derived here imply that we can not avoid having a 
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maximal atmospheric neutrino mixing angle 6 by using NSI with non-zero (p, despite the 
fact that the value of is essentially unrestricted. The reason for this lies in the fact that 
the allowed magnitude of neutrino NSI measured by F is so constrained (due to the lack of 
energy dependence of the NSI evolution equation) that its contribution must be sub-leading. 
This means that a maximum atmospheric neutrino mixing angle is a solid result which must 
be incorporated into any acceptable particle physics model, even in the presence of exotic 
neutrino interactions. 
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APPENDIX: NEW SUPER KAMIOKANDE AND MACRO DATA 

In this section (which does not appear in the published version of this paper) we present 
an update of the hybrid OSC + NSI analysis discussed in Sec. El The calculation of the 
event rates and the statistical analysis is performed according to Ref. [sil, which improves 
the one used so far in essentially three ways: 

Experimental data. Both the Super-Kamiokande and the MACRO collaborations have 
recently released new data. The Super-Kamiokande data used here correspond to 
1489 days (52|, and include the e-like and //-like charged-current data samples of sub- 
and multi-GeV contained events (10 bins in zenith angle), as well as the stopping 
(5 angular bins) and through-going (10 angular bins) up-going muon data events. 
From MACRO we use the through-going muon sample presented in js^, divided in 10 
angular bins. 

Statistical analysis. We now take advantage of the full ten-bin zenith-angle distribution 
for the Super-Kamiokande contained events, rather than the five-bin distribution em- 
ployed previously. Therefore, we have now 65 observables, which we fit in terms of the 
four relevant parameters Am^, 6, e and e'. 
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FIG. 10: Same as Fig. |31 but using the latest Super-K and MACRO data. In the theoretical 
calculation of the expected event numbers the non-zero scattering angle between the incoming 
neutrino and the scattered lepton directions is now properly taken into account. 



Theoretical Monte-Carlo. We improve the method presented in Ref. |12| by properly tak- 
ing into account the scattering angle between the incoming neutrino and the scattered 
lepton directions. This was already the case for Sub-GeV contained events, however 
previously we made the simplifying assumption of full neutrino-lepton coUinearity in 
the calculation of the expected event numbers for the Multi-GeV contained and up- 
going-// data samples. While this approximation is still justified for the stopping and 
thru-going muon samples, in the Multi-GeV sample the theoretically predicted value 
for down-coming i/^ is systematically higher if full collinearity is assumed, as can be 
clearly seen from the second panel of Fig. El The reason for this is that the strong 
suppression observed in these bins cannot be completely ascribed to the oscillation of 
the down-coming neutrinos (which is small due to small travel distance). Because of 
the non-negligible neutrino-lepton scattering angle at these Multi-GeV energies there 
is a sizable contribution from up-going neutrinos (with a higher conversion probability 
due to the longer travel distance) to the down-coming leptons. However, this problem 
becomes less visible when the angular information of Multi-GeV events is included in 
a five angular bins presentation of the data, as previously assumed jlT^ . 

Our results are summarized in Figs. fTUl ITT] and IT^ As already found in Sec. llVl using the 

old data set, the pure NSI solution Am^ = gives a very poor fit, and it is completely ruled 
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FIG. 11: Allowed regions in the e-e' parameter space for the hybrid OSC + NSI mechanism, using 
the new Super-Kamiokande and MACRO data. The oscillation parameters Am^ and sin^(20) are 
integrated out. The shaded areas refer to the 90%, 95%, 99% and 99.73% CL with 2 parameters, 
and the best fit points are indicated by stars. 

out. This occurs since the NSI mechanism is not able to reconcile the anomaly observed in 
the upgoing muon sample with that seen in the contained event sample, as can be clearly 
seen by looking at the NSI line in Figs. El and Hm Conversely, the pure oscillation solution 
e = e' = is in very good agreement with the experimental data. Therefore, we can expect 
that, when combining the two mechanisms of Vy, Vr transition, oscillations will play the 
role of leading mechanism, while NSI's can only be present at a sub-dominant level. 
The global best fit point occurs at the parameter values: 

sin^(2^) = 1, Am^ = 2.3 x 10"^ eV^ e = 6.7 x 10-^ e' = ±1.1 x 10-^ (A.l) 

and it is interesting to note that the new data favour a small but non- vanishing component 
of NSI. From Fig. El we can see that this preference originates from the most vertical 
events in the the Super-Kamiokande and MACRO thru-going /i data samples, which slightly 
favour a stronger suppression of the neutrino signal. However, this effect is not statistically 
significant: the best pure oscillation solution e = e' = 0, which occurs at 6' = 45° and 
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FIG. 12: Behaviour of Ax^ for the new Super-Kamiokande and MACRO data, as a function 
of the oscillation parameters Am^ and 6 (left panels), and of the NSI parameters e and e' (right 
panels). In each panel the undisplayed parameters are integrated out. In the left panels both the 
pure oscillations case and the hybrid OSC + NSI mechanism are displayed. 



/S.m? = 2.5 X 10^'^ eV^, exhibits a which is worse than the global one only by 2.4 
units. The determination of the oscillation parameters Am^ and 6 is very stable under 
the perturbation introduced by the additional NSI mechanism: from the two left panels of 
Fig. ^1 we see that the range of 6 is essentially unaffected, and the only effect of allowing 
for NSI is to slightly weaken the lower bound on Am^. For these parameters we derive the 
ranges 0.84 < sm\2e) < 1 and 1.0 x lO^^ eV^ < Am^ < 4.8 x lO^^ eV^ at 99.73% CL. The 
bounds on the NSI parameters derived in Sec. |3 are strongly improved by the inclusion of 
the new Super-Kamiokande data, and at 99.73% CL we now have —0.03 < £ < 0.02 and 
\e'\ < 0.05. In addition, from the right panels of Fig. Elwe see that the function is quite 
flat in the e' directions for e' 0, and almost symmetric under the exchange e' —e'. 
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